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The two-dimensional Ising model with two kinds of bonds—ferro and antiferromgtgnetic— 
is investigated in both quenched and annealed approximations. High temperature series ex¬ 
pansion technique is used to get the respective free energies and susceptibilities. From the 
analysis of the ratio method results, the phase diagrams are drawn and critical concentrations 
are estimated. It is shown that in both approximations the system behaves similarly, although 
in the annealed case the magnetically ordered phases are more stable against increasing 
impurity concentration. It is also proved that for Bethe lattices initial susceptibilities for 
quenched and annealed systems coincide. 


§ 1. Introduction 

Systems with impurities seem to have attracted the attention of both experi¬ 
mental and theoretical physicists in late years. Because of mathematical difficulties 
mainly the Ising model is considered. The impurities may have magnetic or non¬ 
magnetic (diluted magnetism) character. In either case they may be located at 
the lattice sites (the so-called site problem), or they may change the strength 
and/or sign of the bonds coupling magnetic atoms in the crystal (the bond pro¬ 
blem) . 

As regards the contribution of the impurities to the statistics, two extreme 
cases are considered, the quenched and annealed systems. In the former, impuri¬ 
ties are distributed completely randomly throughout the crystal. They are not 
free to move, but rather “frozen” at their positions. To calculate the free energy 
of a quenched system, one has to compute the free energy for a given distribution 
of impurities, and then average all their possible configurations. In an annealed 
system the impurities (site or bond) are free to move, and their location corresponds 
to the real thermodynamical equilibrium of the system of spins and impurities. 
Therefore, in order to get the free energy of an annealed system, we have to sum 
up simultaneously over spins and impurities configurations. 

The case of diluted magnetism has been studied by many authors (see, e.g., 
Matsubara 1 * and references therein) in both quenched and annealed systems. As 
for magnetic impurities, the annealed approximation (bond system) was mostly 
considered by Kasai, Miyazima and Syozi. 2) ’ 8) The decorated bonds technique was 
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used, which reduced the problem to the case of “pure” Ising model with re¬ 
normalized parameters. 

One of the main questions when dealing with systems with impurities is to 
establish the critical concentration, i.e., the concentration below which there is no 
spontaneous ordering even at zero temperature. In a diluted system there is only 
one such concentration. If we have, however, two kinds of bonds—one ferro and 
the other antiferromagnetic, then there are two critical concentrations, one for 
each type of bond. Rapaport 4) considered extensively the problem of non-magnetic 
impurities in quenched and annealed, as well as in site and; bond approximations. 
Katsura and Matsubara 5) investigated magnetic properties of quenched systems with 
magnetic and non-magnetic impurities* of the bond and site types. 

In this note we examine a two-dimensional Ising model with bond impurities 
(spin 1/2) in both quenched and annealed approximations and find the respective 
phase diagrams. To the best of our knowledge such a comparison, for two-dimen¬ 
sional ferro-antiferromagnetic mixture, has not yet been done. We shall use the 
high temperature series expansion technique and ratio method. 

§2. Free energy of quenched and annealed systems 

The Hamiltonian of the bond problem for an arbitrary lattice has the form 

( 1 ) 

<ij> i 

where is the exchange integral between nearest neighbouring sites i and j, 
and may be either J A (ferromagnetic), or J B (non-magnetic or antiferromagnetic).' 
U is the external magnetic field directed along the z-axis, m — \gg B and (?i= ±1 
is the (reduced) spin operator at the site i. 

We shall now consider more closely the free energies of the quenched and 
annealed systems. Let us denote one of the bond configurations by We label 
the Hamiltonian (1) as when the bond configuration is fixed at $. As 

mentioned in the Introduction, the free energy of the quenched system can be 
written as 

K= 2 \~kTi n Tr «p(-^)} ]n*= £ F($)/N* , (2) 

where Xf denotes the summation over all bond configurations, Tr is the trace 
over spin configurations, F($) is the free energy of a system with a definite bond 
configuration f and N* is the total number of bond configurations 

~ 2 l = (z/2)ffC( z / 2 )jyj, A , (3) 

where z is the number of nearest neighbours, N is the total number of lattice 
sites, Pa(Pb) is the concentration of J A (J B ) bonds. Hereafter we shall use the 
notation <(•••)> for the average over all bond configurations. Therefore (2) may 
be written as 
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F (l =I±F(£)/N* = <F<i£)>. ( 4 ) 


Equation (2) can be expanded, 

F q = — kT <ln{(ch K&) {z/2)Np & (ch X’ B ) (z/2) ^ PB (ch L) N 


X Tr n(l + W/ th Ku) (1 + (Ji th Z/)})>, 

<ij> i 


where 


K i} = ^- (k a = ^-, K b = 

' 2kT\ 2 kT 


2 kT 


, L = ^. 
kT 


( 5 ) 

( 6 ) 


Since we shall derive the phase diagrams from the convergence of the initial 
susceptibility series, we need only second-order terms (thZ,) 2 . We use the ex¬ 
pansion <ln(l+ /)> = </) — -. Taking the trace over spins, expanding 

the logarithm in (5) and truncating the series at (th.L) 2 terms, we get 

F a =F 0 - kT {<E □ > + <£ CD > + <E CZD > + <2 □ □> 

-KQ3Q) (E □)> + ... 

+ OhF) 2 [<E X —x> + <£ x-x>+ ... + <£ Lj-|_j- x > + ... 

+ <(E O x— x»-<(E □) (E x—x)> + ... 

+ < (E CD X-x ) > - < (E CII) (E X-x )> + ...}, (7) 

where 

F 0 =-kT In[2*(ch K A y z ^ Np A ( c h K B )^ Np B( c h L) N ] (8) 

and, e.g., Ll LF* denotes all connected graphs of eighth order in (thi£) 

ending in two “field-points” marked by a cross x, (Xj x - x ) denotes the 

summation of graphs of second order over all their possible configurations and 

(2 Q x — x )« denotes summation of graphs □ and x-x which are not connected. 

On the other hand, the free energy of the annealed system can be written as 

F a = —kTIn |e Tr exp(-^p)} = —kT In Z, (9) 


where Z is the partition function of the equilibrium state. Equation (9) may be 
expanded yielding 

F & = —kT, ln{(eh K A )^ 2 ) Np a( c h K b )^ Np b ( c h L) N 

x S Tr n (! + Wj th K i} ) n (1 + ffi th L) }. (10) 

i <iy> i 

Since 

In 2 (1 +/) = In iV* + In (1 + </», = In iV* + </> - K/> 2 + * * * 
then Eq. (10) can be expressed as 


(ii) 
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F.=F t -*Tln.JV*-*T«SD> + ...+.<J]nn>-i<2]Q><Sa>+'... 

+ (thL) 2 [<sx—x>+...+<2 Ln_r x >+... 

+ <E □ X—x>-<ED)<E X—X> + ... 

+ <U CD x-x>-<X] IZDXE x -x> + (12) 

We notice that F a has an additional term — kTlnN*, as compared with F v which 
is the entropy of bonds arrangement in the annealed system. 


§ 3. Linear chain 

Before we discuss the square lattice, we must mention the distinctive, property 
of the linear chain lattice. Since here are no closed graphs, such as □, we can 
drop all closed graphs in (7) and (12). Consequently, the free energies of the 
linear chain are 


F q =F 0 -^T(thL) 2 {<I] x—x> + <I] x-x> + | 

F a =F 0 -£TlnJV*-£T(thL) 2 K2 x—x> + <I] x-*> + ...}. J (13) 

As the initial susceptibility x is proportional to (th L) 2 in F, both y q and y a for 
the linear chain are identical. The same is also valid for any Bethe lattice. ^ 
for the linear chain is easily calculated from (13): 

ra 2 1+ (p A thK A +p B th Kb) .... 

kT 1 (p A th K a + pQ th K b ) 

This is in agreement with that by Katsura and Matsubara. 6) Terms proportional 
to the fourth- and higher-order powers of (th L) are not identical in and F a , 
even for the chain. 


§ 4. Square lattice 

For the square lattice we may write the free energies of the quenched systems 

F q =F 0 - kTN{Ai + 2 Af+7 Af - 2 AM - Uz* 

+ (th L) 2 [2A X + 64/ + 18 AS + 50 AS + 1424/ + 390 A,* + 1082 A? 

+ 29504/ - 4AM —16AM - 4AM 2 - 7 6AM - 1 6AM* 

- 4 AM 3 ~ 272 AM - 76AM 2 -1 6AM S ] } (15) 

and annealed systems 

Fa=Fo—kT InN* — kTN |4/ + 24/ + ^4/ - 4 AM + (th L) 2 

X [2Ai + 64j 2 + 184/ + 504/ +1424/ + 3944/ +11064/+30584/ - 44/4 2 
- 24 AM ~ 120 AM - 472 AM] J, (16) 


where 
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and 


10A th 2 kT 9 


ICB=th A.. 
2kT 


(17) 


(18) 


Now we can eliminate the Kb bonds from Eqs. (15) and (16), which give F a and 
Fq as functions of the external field H, the ratio x of the coupling constants 

„_Jb 


X- 


Ja 


(19) 


concentrations of the respective bonds p A and p B and the temperature T. Dif¬ 
ferentiating thus obtained F a and F q twice with respect to the external field we 
get the respective susceptibilities and ^ q in the “ferromagnetic-rich” region: 

Xa = T^ ^biKA + 12b? ka + 4 (9b? -f a 3 ) /c A 3 + 4 (25 b? + 6&!<z 3 ) Ka 
kl 


+ 4 (Jib i + 21b?a z - 2b?b 2 + a 5 ) k^+ 4 (195^ 6 + lOOftfo + 6b x a 5 + 3 a? 

- 8b?b 2 -2b?b?) k a 6 + 4 (541 b? + 355b?a 3 + 27 b? a* + 27 b x a? + a 7 
-6b?a z b 2 -2b?b ,- 3 8b?b 2 - 8b?b?~2bJ>?)K a 1 + 4(1475^ 8 + 1170b?a z 
' 4- 150ii 2 a 8 2 + 100^i 3 <2 5 + 6&!<z 7 + 6a z a$ — 32b?a s b 2 — 8b?b± — Ab^ajb? — 4 b?b 2 b 4c 

-13 6b?b 2 - 38b?b? - 8b?b?) a; a 8 } , (20) 


* a = f? {1 + 4il * A + + 4 ( 9 ^ 3 + «•) ^A 3 + 4 (25V + 66,a s ) Ka 4 

+ 4 (71V + 27b i a,, +a t - 26/6 2 ) ff A 5 + 4 (1976/ + 1006/a s + 66,a 6 + 3a s 2 

- 126/6*) k a 6 + 4 (553V+3556/o» + 27 6/a, + 276,0,* + a 7 -66/6*a, 

- 26/6, -606/6,) Ka 7 + 4 (16546/ + 11826/a, + 1906/a, + 1506, V + 66,a 7 

+ 6a s a 5 — 486, 3 6 2 a 3 —126/6, — 2366 / 62 ) Ka 8 } , (21). 


where 


di = x, <2*=i^B«i(l — V), a, = —a 3 (3 — 2a/), 
3 5 


a 7 = —a 6 (45 - 53a, 2 + 17a/), 
105 


bi=pA-\a x Ps, b 2 =p A J r aip B , 6, = 2a,a 3 . (22) 

v e 

Eliminating /C A from £qs. (15) a.nd (16) and introducing the staggered field we 
can get an analogous result for the staggered susceptibilities in the “antiferro- 
magnetic-rich” region. Thus, in both areas ^ has ^ the form of a series in Ka (or 
Kb), and to establish the location of the critical points we can use the ratio method. 
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The phase diagrams deduced from it for the quenched and annealed systems are 
presented in Figs. 1 and 2, respectively, where the critical temperatures are plotted 
versus the concentration p A (=p) of ferromagnetic bonds. In both cases (quenched 
and annealed) we change the valufe of the ratio x by changing the J B coupling 
constant and keeping the ferromagnetic one constant. 



Fig. 1 . Critical temperature versus concentra¬ 
tion p A for the quenched system. The 
right side is ferromagnetic region and the 
left side is antiferromagnetic region. The 
numerical values attached on the curves 
are J B /J A . 
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Fig. 2. The annealed system. The explana¬ 
tion is the same-as that in Fig. 1. 


§5. Discussion 

From inspection of Figs. 1 and 2 it is seen that the effect of impurities is 
in lowering the critical temperature. For the case of equal strengths of ferro 
and antiferromagnetic couplings the curves in both regions are perfectly symmetric. 
With increasing amount of impurities, the temperature of the transition to the 
paramagnetic state decreases, and finally, at a certain concentration, the lattice 
ceases to be magnetically ordered, even.at T=0. This gives the location of the 
critical concentration. It depends strongly on the ratio x. For x = 0 (no antiferro¬ 
magnetic bonds—diluted ferromagnet) the critical concentration for the quenched 
system is (p A =P) A = 0.52 and for the annealed A = 0.47. The rigorous value 
for the annealed system is A = 0.5. 2),s) For all smaller values of p the system is 
in the paramagnetic state. With increasing strength of the antiferromagnetic coupl- 
ing there appears antiferromagnetic region expanding both towards higher critical 
temperatures and higher critical concentrations. On the opposite side, the ferro¬ 
magnetic region “shrinks” towards lower critical concentrations. Since we keep 
J A constant the critical temperature for p = l remains the same. The situation is 
symmetric, i.e., would we" keep J B constant, and use y=x~\ instead of x, we would 
get a picture with ferro and antiferromagnetic curves interchanged. 

From the comparison of the respective curves in Figs. 1 and 2 it is seen that 
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the magnetically-ordered phase is preserved longer in the annealed than in the 
quenched system. The general character of the curves is essentially the same. 
The biggest difference is for the x = 0 curves. It may be concluded that allowing 
the (magnetic or non-magnetic) impurities to move, and therefore to participate 
in establishing the real thermodynamical equilibrium of the lattice, helps in main¬ 
taining magnetically-ordered state. This difference between quenched and annealed 
systems, which vanishes for one-dimensional models, should be even more pronounc¬ 
ed in three-dimensional lattices. 

Since our calculations are high-temperature, series expansion, we do not expect 
the results at low temperatures to be correct. Therefore, we have not obtained 
the plural critical temperatures as they were derived rigorously by Kasai et al. 2),8) 
for the annealed system. However, our results are very accurate at high temper¬ 
atures. 
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